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ABSTRACT 


The  recently  described  phenomenological  model  for  strongly  adiabatic  electron 
transfer  (et)  in  the  lowest  excited  state  Si  of  donor/acceptor  molecules  has  been  developed  in 
more  detail.  The  model  has  been  generalized  to  include  simultaneously  two  types  of  electron 
transfer,  namely  adiabatic  et  in  Si  and  diabatic  et  (S2  -*  Si).  A  general  definition  of  the  et 
"reactant"  survival  probability  is  given  in  terms  of  the  surviving  charge  on  the  donor.  A 
relationship  between  the  survival  probability  and  integrated  emission  intensity  from  Si  and 
S2  is  established.  The  use  of  this  model  is  demonstrated  employing  the  recent  simulation 
results  from  our  group  on  the  excited  state  et  of  bianthryl  (BA)  and 
4-{ 9— an t hryl )-N , N * — di  met hy lanili ne  ( AD M  A ) . 
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1.  Introduction 

Photophysical  studies  on  organic  aromatic  electron  donor/acceptor  compounds  are 
leading  to  new  information  on  several  issues  of  contemporary  interest,  including:  the  role  cf 
solvation  dynamics  in  electron  transfer  (et)  reactions;  the  mechanism  of  ultrafast  charge 
separation  (which  is  relevant  to  problems  ranging  from  photosynthesis  to  nonlinear  optics) 
and  the  effect  of  the  solvent  coordinate  on  the  absorption  and  emission  spectroscopy  of  the 
aromatic  compounds  in  polar  solvents.1  The  comparison  between  experiment  and  et  theory 
has  been  a  particular  emphasis.  Historically,  one  of  the  most  extensively  studied  organic 

2 

donor/acceptor  compounds  is  4— (9— anthryl)— N,N’— dimethylaniline  (ADMA)  (Figure  1).“ 
The  preceding  article  in  this  issue  presents  static  and  ultrafast  spectroscopic  data  on 

n 

ADMA.  The  results  are  analyzed  in  terms  of  a  phenomenological  et  model  (in  the  strongly 
adiabatic  regime)  that  was  recently  introduced  by  Kang  et  al.  for  et  between  strongly 
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coupled  donors  and  acceptors. 

In  this  paper  we  explore  the  strongly  adiabatic  theories  in  more  detail.  The  strongly 
adiabatic  approach  is  apparently  a  more  accurate  description  of  et  between  strongly 
coupled  donors  and  acceptors  than  the  more  conventional  et  theories  which  are  appropriate 
for  weak  to  moderate  electronic  coupling.1’^  111  It  is  important  to  note,  however,  that  all  of 
these  theories  invoke  the  diabatic  (localized)  states  of  Marcus  theory,  here  denoted  fA  and 
♦  B,  where  A  and  B  refer  to  the  reactant  and  product  states,  respectively.1^  The  magnitude 
of  the  interaction  between  the  localized  states  is  given  by  the  electronic  coupling  matrix 
element,  Hin  =  f  t  /ftndr.  The  et  kinetics  are  induced  by  fluctuations  and  vibrations  of 
the  nuclear  coordinates,  including  the  solvent  and  intramolecular  vibrational  modes.  In  the 
nonadiabatic  and  moderate  adiabatic  theories  H.D  is  treated  by  perturbation  theory.  In  the 
presence  of  a  large  HAD  the  perturbative  treatment  is  inappropriate  and  a  strongly  adiabatic 
model  becomes  more  useful.  A  rigorous  theory  of  the  coupling  between  the  solvent 
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coordinate  and  the  electronic  coordinate  in  the  strongly  adiabatic  regime  has  been  given  by 

Kim  and  Hynes.  ^  Molecular  dynamics  simulation  models  that  are  compatible  with  the 

12 

strongly  adiabatic  regime  have  also  recently  been  published. 

3 

For  the  specific  case  of  excited  state  et,  Tominaga  et  al.  showed  the  electronic 

coupling  between  the  locally  excited  (LE)  state  and  charge  transfer  (CT)  state  of  ADMA 

13 

(Figure  1)  is  so  great  that  it  is  more  useful  to  use  the  approximate  adiabatic  states  Si  and 
S2  rather  than  the  diabatic  states.  Figure  2  portrays  the  energy  dependence  of  the  diabatic 
states  on  the  displacement  of  the  so-called  solvent  coordinate,  z.  Figure  2  shows  the  energy 
dependence  of  the  adiabatic  states  which  are  mixtures  of  the  diabatic  states. 

|S,>  =  C<  ”(z)|LE>  +  C“>(z)|CT>.  (1-1) 

|S2>  =  C<  2)(Z)|LE>  +  C<2)(z)|CT>.  (1-2) 

The  various  parameters  required  to  construct  Figure  2  were  determined  empirically  by  a 
phenomenological  approach  which  involves  fitting  the  experimental  spectra  to  theoretical 

4 

predications  using  the  curves  in  Figure  2  and  other  information,  see  below  and  Kang  et  al. 
for  more  details. 

For  ADMA  the  diabatic  states  involved  in  the  excited  state  et  are  the  LE  state  which 
*  * 

results  from  a  t-t  transition  of  the  anthracene  ring  (DA  )  and  the  CT  state  which  is  an 
intramolecular  ion— pair  state  (D^A- )  (Figure  1). 

The  ultrafast  spectroscopy  of  ADMA  shows  several  features  which  are  characteristic 
of  an  adiabatic  mechanism,  and  would  be  awkward  to  treat  in  a  nonadiabatic  representation. 
In  particular,  the  fluorescence  spectrum  is  structureless  and  exhibits  a  continuous  shift  of 
emission  maximum  frequency  on  a  time  scale  which  is  approximately  equal  to  the  known 
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time  scale  for  solvation  dynamics,  e.g.  1.7  ps  in  N.N-dimethylformamide.  ’  The  decay  of 
the  emission  intensity  is  not  due  to  loss  of  population  of  excited  molecules,  which  occurs  on  a 
much  slower  time  scale.  The  observed  behavior  is  close  to  simulations  using  the  adiabatic 
model.  In  terms  of  this  model,  the  fluorescence  is  assigned  to  the  Si  state.  The  drop  in 
intensity  and  evolution  of  the  emission  maximum  frequency  is  assigned  to  diffusion  along  the 
solvation  coordinate  which  causes  the  probability  distribution  of  the  solvent  to  evolve  from 
its  initial  photoprepared  distribution  (peaked  near  z  =  0)  to  its  equilibrium  distribution  in  Si 
(peaked  near  z  =  1). 

The  decrease  in  intensity  as  the  solvent  coordinate  evolves  is  due  to  the  adiabatic  et 
process.  As  the  average  z  values  goes  from  zero  to  one  the  amount  of  LE  character  in  ♦ 

^  1 

decreases,  and  the  amount  of  CT  character  in  f0  increases.  The  transition  moment  (St  - 
S0)  simultaneously  weakens,  which  causes  the  integrated  fluorescence  intensity  to  decrease. 

The  Si  adiabatic  et  processes  of  ADM  A  and  related  examples  are  beyond  the 
conventional  definitions  of  et  kinetics.  It  is  not  immediately  obvious  how  the  rate  coefficient 
(ket)  for  this  process  should  be  defined  in  theoretical  or  experimental  terms.  This  general 
problem  is  the  main  focus  of  the  present  paper.  In  addition,  the  underlying  physical 
principals  and  approximations  of  the  phenomenological  adiabatic  excited  state  et  model  are 
delineated.  An  additional  focus  of  this  paper  is  the  existence  of  two  distinct  et  processes,  i.e. 
the  S2  Si  "nonadiabatic"  et  mechanism  and  the  Sj  (or  S2)  adiabatic  et  mechanism,  which 
were  not  classified  in  the  original  Kang  et  al.^  paper.  This  paper  is  based  primarily  on 
simulations  within  the  adiabatic  model  on  ADMA  and  BA.  Some  of  the  simulation  results 
have  been  previously  discussed  and  shown  to  be  in  semi  quantitative  agreement  with 
experiment.  In  the  present  paper,  we  show  that  the  simulation  results  can  be  analyzed  via  a 
general  definition  of  the  rate  coefficient  to  extract  kinetic  information  on  the  et  process. 
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2.  The  Adiabatic  Electron  Transfer  Model  and  ket 

This  section  develops  the  basic  adiabatic  et  model  with  particular  emphasis  on  the 
various  inherent  approximations.  We  restrict  this  discussion  to  molecules  like  ADMA  which 
involve  a  "two  state"  et  process,  i.e.  LE  -♦  CT  in  a  diabatic  picture.  The  molecule  BA  is 
more  complicated  since  there  are  two  CT  states  and  LE  stated’**  However,  the  theoretical 
results  for  ADMA  are  applicable  to  BA  with  a  simple  extension,  as  will  be  discussed  in 
Section  5. 

For  purposes  of  describing  the  excited  state  et  of  ADMA  and  the  photodynamics  of 
this  compound,  it  is  important  to  consider  the  So  ground  electronic  state.  While  this  state 
plays  a  small  direct  role  in  the  LE  ->  CT  et,  it  does  play  an  essential  role  in  the  overall 
photodynamics.  Optical  excitation  from  So  to  Sj  and  S2  prepares  the  t  =  0  configuration  of 
ADMA,  and  transient  emission  from  the  excited  ADMA  to  So  is  an  observable  which  can  be 
used  to  measure  ket  and  other  dynamical  observables  (see  below).  For  the  purpose  of  this 
paper  we  will  assume  that  LE  and  So  have  identical  charge  distributions,  which  is  consistent 
with  experiment,  ’  such  that  both  states  have  identical  equilibrium  solvent  configurations. 

2.1.  Diabatic  States  versus  Adiabatic  States 

Invoking  Marcus  theory, Tominaga  et  al.^  have  shown  that  the  free  energy 
dependence  of  So,  LE,  and  CT  should  be  given  as  follows: 


F„(»)  =  (l/2)ksz2, 

(2-1) 

Fle(z)  =  (l/2)kszJ  +  Fo, 

(2-2) 

FCTM  =  (l/2)ks(Z— l)’+Fo  — 4F«, 

(2-3) 
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where  ks  is  the  solvent  force  constant,  F0  is  the  energy  of  the  spectroscopic  transition 
between  So  and  LE,  and  AF°  is  the  equilibrium  free  energy  change  for  the  LE  -*  CT  charge 
separation. 

It  is  useful  to  consider  approximate  stationary  adiabatic  states  for  the  ADMA 
problem  as  represented  in  Figure  2  and  by  eqns  1-1  and  1-2.  The  z  dependent  energies  of 
these  two  states  can  be  found  approximately  by  diagonalizing  the  secular  equation, 


flE(z>-f 

U 

LE,CT 


TJ 

LE,CT 

fctW-f 


(2-1) 


where  the  diagonalized  free  energies  are  given  by  eqns  2—2  and  2—3,  H  is  the  electronic 
coupling  between  LE  and  CT,  and  F  is  the  adiabatic  free  energy.  Accordingly,  the 
coefficients  can  be  found  by  the  usual  means. 


2.2.  Electron  Transfer  Kinetics  and  the  Nonadiabatic  and  Weak  Adiabatic  Theories 

The  basic  formulation  of  et  theory  in  solution  is  firmly  routed  in  the  localized  two 
state  description  of  Marcus,  which  involves  a  reactant  state  fR  and  product  state  f  p.10  In 
the  language  of  this  paper  the  electron  that  is  being  transferrer  in  the  reaction  is  first  on  the 
donor  group,  localized  in  fLE,  and  is  later  on  the  acceptor  group,  localized  in  fCT-  As  in 
Huckel  theory,  the  overlap  between  the  localized  states  is  negligible. 

<LE|CT>  =  0.  (2-5) 


Based  on  this  picture  one  can  envision  an  operator  that  measures  the  charge  on  the  donor 
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such  that 


<LE  | 

«D  1  LE>  2  qei, 

(2-6) 

<CT 

l?c|CT>  2  0. 

(2-7) 

Thus,  transitions  between  LE  and  CT  physically  transfer  charge  between  the  donor  and 

acceptor.  Many  of  the  variations  of  contemporary  et  theory  are  based  on  an  approximate 

treatment  of  the  coupling  of  the  nonadiabatic  states  to  the  medium  fluctuation,  via  a 

Liouville  equation  description,  that  involves  the  diabatic  states  (LE  and  CT)  as  a  basis  set 

5—10 

and  a  two  state  Hamiltonian: 


Hei(f)  = 


f  Hle,ct 

K.E,ct1E<1>  J 


(2-8) 


where  t  is  the  energy  "exchanged"  with  the  quantum  mechanical  nuclear  and  electronic 
degrees  of  freedom,  and  AE(t)  is  a  fluctuating  energy  that  describes  the  coupling  and  thermal 
motion  of  the  classical  degrees  of  freedom. 

The  Liouville  equation  involves  a  two  by  two  density  matrix  7  in  the  diabatic  basis 
set,  |  LE>  and  | CT>: 

3^  =  £[iye).7(</t)l.-  (KlT(f/t)]..  (2-9) 


Here  we  employ  the  notation  [A,B]±  =  AB  ±  BA.  K  is  a  stochastic  relaxation  operator 
which  induces  relaxation  in  the  product  (CT)  state. 

The  solution  to  this  equation  typically  involves  several  assumptions:  (i)  The  initial 
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fl  Ql 

concentration  is  ?(f/0)  =  g  g  =  |LE>.  (ii)  AE(t)  is  the  fluctuating  energy  gap,  i.e.  AE(t) 
=  FL£(z(t))  —  FCT(z(t))  (see  eqns  2—2  and  2—3),  where  z(t)  is  often  assumed  to  obey  a 
Fokker— Planck  equation  for  stochastic  (Brownian)  solvent  coordinate  motion  on  harmonic 
potential.  (Although  we  have  emphasized  only  the  solvent  coordinate  in  the  reaction 
coordinate,  the  theory  can  be  extended  to  include  the  intramolecular  vibrational  motion.) 

Several  of  the  recent  theoretical  results  on  et  kinetics  that  include  the  effect  of 
solvation  dynamics  (and  in  some  cases  vibrational  modes)  are  closely  related  to  the 
perturbative  solutions  of  eqn  2— 9.^  ^  Since  the  physical  identity  of  the  |  LE>  and  jCT> 
states  is  maintained  in  this  nonadiabatic  approach,  the  reactant  survival  probability  is 
simply  given  for  a  solvent-coordinate-only  model  by 

SNA(0  =  /  dzpLE(z,t),  (2-10) 

where  pL£( z,t)  is  the  time-dependent  classical  probability  distribution  of  the  reactant  in  the 
solvent  coordinate.  The  reaction  rate  coefficient  is  conventionally 

ket(t)  =  -dlnSNA(t)/dt.  (2-11) 

Note  that  SyA(t)  is  simply  proportional  to  the  surviving  electronic  charge  on  the  donor  (see 
eqn  2-6). 

2.3.  Electron  Transfer  Kinetics  in  the  Strongly  Adiabatic  Regime 

For  the  excited  state  et  reactions  involving  a  donor/acceptor  pair  linked  by  a  single 
bond,  the  electronic  matrix  element,  H  is  large  enough  that  the  perturbative  treatment 

LLjL  i 

in  the  nonadiabatic  and  weak  adiabatic  regime  becomes  questionable.  In  particular,  the 


degree  of  adiabaticity  can  be  quantified  by  the  adiabaticity  parameter, 


H,= 


8tH  3r 
LE,CT  L 


ftk, 


(2-12) 
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where  r  is  the  longitudinal  dielectric  relation  time.  For  ADMA  and  BA  in  common  polar 
solvents,  Ha  »  1  (~  80  with  ks  =  20  kcalmol'1,  HL£  CT  =  1  kcalmol'1,  and  TL  =  1  ps), 
indicating  that  strong  adiabatic  effects  should  be  considered.  In  addition,  the  adiabatic  free 
energy  surfaces  for  ADMA  have  no  barrier  separating  the  initially  prepared  (reactant-like) 
and  equilibrium  (product— like)  regions.  This  is  a  strongly  adiabatic  effect  which  is 
inconsistent  with  a  perturbative  treatment  of  H  .  Hynes  has  considered  strong 

LLjL  i 
6 

adiabatic  et  reactions  over  a  substantial  barrier,  and  related  this  regime  to 

15  19 

transition-state— theory.  But,  our  concern  here  is  small  barrier  and  barrier— free  ’ 
dynamics  on  the  adiabatic  surfaces. 

We  have  been  particularly  interested  in  an  "outer-sphere"  strongly  adiabatic  et 
model,  which  assumes  that  vibrational  coordinates  play  a  minor  role,  and  treats  the  solvent 
coordinate  z  classically  at  each  value  of  z.  The  system  can  be  represented  by  a  two  by  two 
density  matrix  f(z,t), 


r(z,t)  =  />( z,t) » 

q>> W 1  +cy> ( i ) !  q>> (»)c <  j ■  (z)+c< , (*) 

cU'McLl 1  (^)+c' ■' (z)qz, (z)  q>> (z) 1  +c< >> ( . ) ! 


(2-13) 
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Here  the  off-diagonal  elements  are  geueiaily  time— dependent  and  signify  coherence  between 
the  Si  and  S2  states.  p(z,t)  is  the  classical  probability  of  the  solvent  coordinate  at  z. 

It  is  important  to  consider  the  properties  of  the  off-diagonal  elements.  Roughly 
speaking,  motion  along  z  will  produce  coherence,  i.e.  non-zero  off-diagonal  elements.  Also, 
the  laser  excitation  pulse  which  prepares  Si  and  S2  from  S0  will  create  "coherence"  among 
So,  Si  and  S2.  However,  we  make  the  approximation  that  the  off-diagonal  terms  in  the 
adiabatic  state  representation  are  2  Chat  all  relevant  times,  (the  random  phase 
approximation).  This  is  appropriate  where  the  S 1/S2  energy  gap  is  large  enough  to  ensure 
rapid  dephasing  compared  with  the  time  scale  for  z  fluctuations.  Under  these  circumstances 
eqn  2-14  applies,  if  it  is  further  assumed  that  motion  along  the  solvent  coordinate  z  satisfies 

1C  1C 

a  generalized  Smoluchowski  equation  (GSE).  ’  In  an  earlier  paper  we  used  a  more  exact 

approach  involving  the  generalized  Langevin  equation  (GLE)  which  is  computationally  more 
4 

intensive.  A  comparison  of  the  GLE  and  GSE  treatments  is  made  in  Section  5. 

+  (««) 

Here  D(t)  and  F(z)  are  a  generalized  diffusion  coefficient  and  a  potential  acting  on  the 
system,  respectively. 


dp  (z,t)  o  o  ,  dT  (z) 

- =  +  HT  W^p  s/z,t)  “  ksrs2(z)/,Si(z,t)  +  ^s/2^2’^' 

(2-15) 

dp  (*,t)  a  a  1  dFc  (z) 

~W  =  +  KTTz-2  +  kS^ZKM  ~  kS2-s,(zK2(z,t)- 

(2-16) 
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Here  kg  (z)  and  kg  (z)  are  defined  as  diabatic  rate  constants  and  represent  "internal 

conversion"  between  S!  and  S:-  In  this  paper  we  leave  the  form  of  these  internal  conversion 

mechanisms  unspecified.  They  are  physically  analogous  to  the  radiationless  decay  between 

states  of  the  same  spin  multiplicity,  and,  as  such,  formally  represent  a 

non— Born— Oppenheimer  transition  of  state.  In  order  to  simulate  the  overall  dynamics  in 

terms  of  eqns  2-15  and  2-16,  initial  distributions  of  p  (z,t)  and  p  (z,t)  are  required,  and 

values  for  the  various  parameters,  D(t),  ks,  ...,  etc.  are  needed  The  methods  for  performing 

such  simulations  from  phenomenological  estimates  have  been  outlined  previously  by  Kang  et 
4  3 

al.  and  Tominaga  et  al.,  and  will  also  be  discussed  below. 

Our  major  interest  here,  however,  is  a  "reactant"  survival  probability  and,  in  turn,  a 
definition  for  ket  that  is  appropriate  in  all  et  regimes  including  the  strongly  adiabatic  limit. 
The  definition  of  the  survival  probability  follows  from  the  definition  of  the  operator  q  (eqns 
2-6  and  2-7).  The  expectation  value  of  this  operator  for  the  system  wavefunction 
(according  to  the  random  phase  approximation)  can  be  shown  to  be 

f  dz<fsyS|qD|tSys>  =  qelt  /  (z)!dz  +  /pSj(  «.t)Cj>2)(«)Sdl). 

(2-17) 

Thus,  the  "generalized"  survival  S„(t)  of  the  electronic  charge  on  the  donor  divided  by  the 
charge  of  an  electron  q£i  is 

sg(‘)  =  /d^SlMcL i'W  +  »s,(*.‘)c<” «’] 

=  (CM>I(t)  +  Cjj>’!(t))-NT  (2-18), 
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where  the  rhs  of  eqn  2-18  represents  the  z-averaged  probability  of  the  LE  character  times 
the  total  number  N_  of  Si  and  S2  molecules.  Nc  and  Nc  ,  respectively. 

T  5>t  i>2 


=  N  +  N  . 

Si  S2 

(2-19) 

=  /  d»PSl(«.‘). 

(2-20) 

,  =  /d^s, («•*)• 

(2-21) 

An  average  of  a  function  of  z,  C(  1J  (z),  over  z  on  the  Si  surface  at  time  t  is  defined  as 

X 

ci4,(t)  =  [  f  dzpsXz, t)C^  (z)]/[f  dzpSi(z,t)+  /dz/?S2(z,t)], 

(2-22) 

where  i=l  or  2  and  x  =  LE  or  CT.  In  turn,  the  generalized  et  rate  coefficient,  ket(t),  can  be 
expressed  as 


ket(t)  =  — dln(SG(t)  -  S0(.))/dt.  (2-23) 

This  expression  is  analogous  to  the  conventional  definition  of  a  reaction  rate  constant  in 
terms  of  a  surviving  "concentrations"  of  the  reactants.  SG(®)  is  the  equilibrium  value  of  the 
survival  of  the  electronic  charge.  In  the  case  of  the  adiabatic  model,  SG(®)  is  the  thermally 
average  LE  character  at  equilibrium.  By  analogy  to  conventional  reactions,  if  SG(®)  is 
approximately  zero,  the  adiabatic  et  reaction  is  effectively  irreversible.  In  contrast,  a 
non— zero  SG(®)  signifies  a  reversible  et  reaction,  which  is  the  case  for  the  Si  et  reaction  of 
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ADMA  (see  below). 

Having  defined  a  general  survival,  we  now  turn  to  its  relationship  to  experimental 
observables.  The  key  observable  is  the  integrated  emission  intensity 

A(t)  at  J di/I(t,i/)/z/3,  (2-24) 

where  I(t ,z/)  is  the  time— dependent  fluorescence  spectrum  to  So  from  both  Si  and  S2  (within 
the  random  phase  approximation).  An  expression  for  A(t)  in  terms  of  SG(t)  can  be  derived 
by  simple  manipulations  of  the  basic  photodynamic  equations  of  the  Kang  et  al.  model.  A(t) 
is  proportional  to  the  square  of  the  expectation  value  of  the  electronic  transition  dipole 
moment  operator  n  (within  the  Franck— Condon  approximation). 

A(t)  a  f  dz |  <fsys|/<|So>  j 2  (2-25) 

«  / dz[/JSi(z,t)|<Si|/i|S0>|2  +  ,S2(z,t)|<S2MSo>|2]. 

(2-26) 

The  z  dependent  contribution  to  the  transition  moment  is  given  by 

<SiU|So>  =  C«  +  C<  ”  (z)„CT,  (2-27) 

where  i  =  1  or  2  and  fi  and  are  the  dipole  transition  moments  connecting  the  diabatic 

L  E  CT 

states, 

Mle  =  <LE|/i|S0>,  (2-28) 
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/*CT  =  <CT|/j|S0>. 


(2-29) 


Utilizing  the  definition  in  eqn  2—22  of  the  average  of  over  the  solvent  coordinate  at  time 
t,  it  is  useful  to  make  the  approximation  that  the  average  of  the  square  is  approximately 
equal  to  the  square  of  the  average,  i.e. 


C<i>2(t)~[C<*>(t)] 


LE 


(2-30) 


and  similarly 


<2-31> 

In  a  latter  section  of  this  paper  we  examine  the  range  of  validity  of  this  approximation  using 
the  simulation  results. 

The  integrated  emission  intensity  A(t)  can  be  related  to  C<^>  and  C^>  by  the 
following  expression  from  eqns  2—26  and  2-27 


AW°(Cti‘hC  +  +  <C1 1'^LE  + 

(2-32) 

If  all  the  population  on  S2  relaxes  to  the  Si  surface  due  to  the  fast  internal  conversion,  A(t) 
can  be  related  to  SQ(t)  in  a  simple  way,  using  eqns  2-18,  2-30,  2-31,2-32  and  (C£  (t))2  + 

(c<‘>(t))2  =  i, 
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(2-33) 


where 

B(t)  =  +  C^M„CI)’,  (2-34) 

A(od) 

where  we  used  NT  =  1  (slow  population  depletion  process  compared  with  the  et  reaction). 

This  expression  is  a  general  definition  for  the  survival  in  all  regimes  from  the  weakly 
adiabatic  to  the  strongly  adiabatic,  within  the  previously  noted  approximations. 

The  situation  is  particularly  simple  in  the  non-adiabatic  limit  for  which  the  weak 
coupling  between  the  donor  and  acceptor  ensures  that  pLE  »  In  this  limit,  the  optical 
excitation  pulse  prepares  a  pure  LE  state,  and  only  LE  -*  So  emission  can  be  observed. 
Under  these  circumstances  eqn  2—33  reduces  to  eqn  2—10. 

SG  (0  *  AWlMo  *  S„»«'  <2-35> 

which  means  that  the  integrated  emission  intensity  directly  monitors  the  LE  concentration. 
Furthermore,  in  the  nonadiabatic  regime  if  the  reactant  and  product  wells  are  divided  by  a 
large  activation  barrier,  the  survival  can  be  simply  measured  by  monitoring  the  emission 
intensity  at  a  single  wavelength  region 

SG(t)  «  I(t,4  (2-36) 
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For  the  strongly  adiabatic  limit  eqns  2—35  and  2—36  are  not  generally  valid.  Indeed 
they  very  poorly  describe  the  situation  for  ADMA.  On  the  other  hand,  eqns  2—35  and  2-36 
can  be  approximately  valid  in  the  strongly  adiabatic  limit  if  a  large  activation  barrier  is 
present.  Under  this  situation,  the  shape  of  /?(z,t)  in  the  reactant  "well"  can  be 
approximately  time— invariant.  These  issues  should  become  clearer  as  we  consider  specific 
simulation  of  the  excited  state  et  dynamics  of  ADMA  and  BA. 

3.  An  Approximate  Expression  of  ket(t)  in  the  Strongly  Inverted  Regime 

Ultrafast  laser  studies  have  been  made  for  BA,  ADMA,  and  related  compounds.  In 
the  case  of  BA  the  excitation  wavelength  in  these  studies  exclusively  prepares  the  Si  state. 
For  ADMA,  both  Si  and  S2  are  prepared  by  the  laser  excitation,  but  the  S2  state  is 
converted  to  Si  within  <  150  fs.  Thus,  for  both  compounds,  the  et  process  occurs 
predominantly  in  the  St  state.  For  et  in  the  Si  state  only,  eqn  2-18  reduces  to 

sG(t)  =  /di(iSi(I,t)q  »(*)*.  (3-1) 

The  functional  form  of  CT(  *>  (z)  is  obtained  by  solving  the  secular  equation  (eqn  2—4). 

L  L 

For  example, 


LE,CT 


t4I*LE,CT 


+  AF(z)J  +  AF(i)J+4HlIjCI2 1", 


(3-2) 


where  AF(z)  is  free  energy  difference  between  the  LE  and  CT  states,  and 
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IF(z)  =  Fle(z)  -  Fct(z) 


=  ks(z  -  1/2)  +  4F°.  (3-3) 

In  general,  SQ(t)  and  correspondingly  ket(t)  can  be  calculated  by  first  employing  the 
GSE  and  GLE  to  calculate  p( z,t)  from  a  given  initial  condition  p(z,t=0).  The  procedure  is 

O  A 

outlined  elsewhere. 

In  order  to  establish  a  basic  physical  understanding  of  SG(t)  and  ket(t)  it  is  useful  to 

consider  an  approximate  analytical  expression  for  these  questions  in  the  limit  of  a  strongly 

inverted  et  reaction,  such  as  ADMA.  In  this  limit,  F  (z)  has  a  single  minimum  at  z  =  1, 

s  i 

and  is  approximately  equal  to  F  (z)  near  the  z  =  1  minimum,  i.e.  F  (z)  is  harmonic  near  z 
=  1  as  shown  in  Figures  2  and  3.  It  will  be  shown  below  that  at  late  times  p( z,t)  is 
approximately  Gaussian  with  a  maximum  position  z  ~  0.95,  that  is  approximately  the 
equilibrium  value,  z  =  1,  see  the  bottom  panel  in  Figure  3.  Under  these  circumstances  a 
solution  to  eqn  2-14  for  a  harmonic  potential  leads  to  the  following  expression  for  z, 

z(t)  2  1— exp(— t/r8),  (3—4) 

where  it  has  been  assumed  that  the  solvation  dynamics  are  overdamped  and  well 
characterized  by  a  single  relaxation  time  rs.  Eqn  3—4  is  the  same  result  found  in  the 
derivation  of  the  so-called  transient  Stokes  shift. * 

Assuming  that  the  shape  of  p( z,t)  is  slowly  varying  on  the  time  scale  of  rs,  which  is  a 
result  of  the  simulation  (Figure  3),  we  can  replace  the  integral  in  eqn  2—18  by  the  following 
expression, 
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ss(0  =  c<»(i(t))J, 


(3-6) 


which  should  apply  at  long  times. 

Combining  eqns  3—2,  3—4,  3—5,  and  the  definition  of  ket(t)  in  terms  of  S„(t)  (eqn 
2-23),  it  is  easy  to  show  that  the  limiting  rate  coefficient, 

kljra  =  Urn  ket(t) 

=  !i,m(-^ln(Se(t)-S0(.))],  (3-0) 

leads  to 

kJi"  =  1/Ts.  (3-7) 

This  relationship  indicates  that  for  adiabatic  et  reaction  in  the  strongly  inverted 
regime,  solvation  dynamics  are  rate  limiting,  which  is  closely  related  to  analogous  results 
from  the  weakly  adiabatic  model.  *’^~ *** 

4.  The  Si  Adiabatic  ET  Kinetics  of  ADMA 

Q 

Tominaga  et  al.  have  reported  transient  emission  data  I(i/,t)  on  ADMA  in  the  polar 

aprotic  solvent  N,N-dimethylformamide  (DMF)  which  according  to  transient  Stokes  shift 

measurements*'*5  has  an  average  solvation  time  <rs>  of  1.0  —  1.5  ps.  The  experimental 

I(t,i/)  results  are  very  similar  to  simulations  of  the  spectra  using  strongly  adiabatic  model,  as 

3  4 

shown  in  Figure  4.  The  simulations  in  Figure  4  have  been  described  previously.  ’  They 
involve  a  convolution  of  an  empirical  spectral  shape  function  with  p(z,t), 
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«  f  dz|  <Si|/i|  So>  |  i/0(z))/)(z,t)j/3, 

(4-1) 

where  g(i/o(z)>I/~I/o(z))  the  normalized  emission  spectral  shape  and  j/0(z )  corresponds  to  the 
free  energy  difference  between  the  St  and  So  states,  v^z)  =  (F  (z)  -  F  (z))/h.  Here  p{ z,t) 

b  1  5>0 

3 

was  calculated  numerically  from  eqn  2—14  using  an  initial,  spectroscopically  prepared 
distribution,  see  the  bottom  panel  in  Figure  3.  In  the  simulation  we  used  a  two  exponential 
function  for  the  solvation  dynamics  C(t),11) 

C(t)  =  Arexp(— t/rj)  +  A2-exp(-t/r2).  (4-2) 

These  values  for  C(t)  were  measured  by  the  transient  Stokes-shift  method,  as  discussed 
elsewhere.  The  time-dependent  diffusion  coefficient  which  is  used  in  the  GSE 
simulations  can  be  related  to  D(t).  For  C(t)  of  DMF  Ai  =  0.55,  A2  =  0.45,  tx  =  0.75  ps, 
and  r2  =  2.5  ps.1  ^  Unfortunately,  the  experimental  measurements  of  the  transient 
fluorescence  spectrum  are  not  reliable  over  a  broad  enough  wavelength  region  to  allow  for  an 
accurate  determination  of  A(t),  the  key  observable  for  the  determination  of  SG(t).  We  are  in 
the  process  of  improving  the  apparatus  employed  in  these  measurements  in  order  to 
determine  A(t)  accurately.  In  this  paper,  however,  we  will  emphasize  the  simulated  A(t) 
data  as  shown  in  Figure  4.  The  simulated  results  (eqn  4—1)  will  allow  us  to  demonstrate  the 
physical  significance  of  the  SG(t)  definition  (eqn  2-33)  and  allow  us  to  evaluate  one  of  the 
important  approximation  (eqns  2-30  and  2-31). 

Tominaga  et  al.  have  demonstrated  that  the  emission  dynamics  of  ADMA  at  times 
>  150  fs  is  primarily  due  to  emission  from  Si.  The  time-dependent  spectra  (Figure  4) 
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reflect  the  continuously  evolving  probability  distribution,  pSi(z,t),  in  St  as  shown  in  Figure 

3.  Note  that  as  the  distribution  shifts  toward  equilibrium  (near  z  =  1)  the  emission 

intensity  shifts  toward  lower  frequency  since  the  So/Si  gap  decreases  as  z  goes  from  zero  to 

one  (Figure  3).  In  addition  to  shift  in  the  spectrum  as  Pg^z.t)  evolves,  the  integrated  area 

(Figure  5)  also  shows  time-dependence.  Simply  speaking,  as  the  average  value  of  z  increases 

from  0  — *  1  the  LE  character  of  Si  decreases,  and,  as  a  result,  A(t)  decreases  because  > 

a  Note  that  the  actual  number  of  Si  molecules  is  essentially  constant  during  the 
^CT 

evolution.'* 

Note  also  that  S_(t)  is  not  zero  at  infinite  time,  since  the  relaxed  form  of  Si  is  not 
purely  CT.  Therefore,  S„(»)  has  a  finite  value  in  eqn  2-23.  ket(t)  is  portrayed  in  Figure  5. 
The  time-dependence  of  ket(t)  reflects  three  complex  effects,  namely  the  dependence  of 
C<  11  (z)  on  z,  the  anharmonic  shape  of  F  (z)  at  the  initially  prepared  z  values,  and  complex 
diffusion  effects  due  to  the  time-dependent  diffusion  coefficient  D(t),  as  shown  in  Figure  4 
Nevertheless,  a  unique  value  for  the  et  rate  constant  can  be  obtained  by  determining  the 
limiting  rate  coefficient  ketm- 

The  generalized  survival  during  the  Si  evolution  is  shown  in  Figure  5.  This  shows 
that  the  survival  probability  indeed  decreases  as  the  Si  adiabatic  et  process  occurs.  The 
value  we  extract  for  ADMA  from  the  simulation  is  kit"1  =  1/(2. 50  ps).  kltm  is  equal  to  the 
inverse  of  the  slower  solvation  time  r 2.  This  is  precisely  the  result  expected  in  the  limit 
discussed  in  Section  3,  where  the  potential  is  approximately  given  F^z)  in  the  equilibrium 
region  of  Si  and  the  solvation  dynamics  are  characterized  by  C(t)  of  the  form  of  eqn  4-2. 

Figure  6  provides  a  comparison  of  the  survival  Sg(t)  and  hence  ket(t)  calculated  two 
ways.  The  middle  panel  shows  ket(t)  calculated  "exactly”  from  eqn  2-23.  Alternatively,  we 
can  determine  S_(t)  from  the  simulated  area  using  the  approximate  result  which  connects 
experiment  and  S  (t)  (eqn  2-26).  The  upper  panel  shows  ket(t)  resulting  from  a  simulation 
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of  S„(t).  For  the  both  cases  a  two  exponential  form  for  the  solvation  was  used.  While  the 
exact  and  approximate  results  are  not  identical,  the  approximate  approach  is  reasonably 
accurate,  especially  at  longer  times,  supporting  the  validity  of  the  approximation  of  eqns 
2-30  and  2-31. 

It  should  be  strongly  emphasized  that  while  ket(t)  can  be  simply  related  to  A(t), 
there  is  no  simple  relationship  between  ket(t)  and  the  emission  dynamics  I(i/,t)  at  any  single 
emission  wavelength.  This  is  an  important  result  which  impacts  on  the  interpretation  of 
published  data  on  the  transient  fluorescence  data  of  ADMA  and  related  molecules. 

The  bottom  panel  in  Figure  6  shows  "exact"  ket(t)  (eqns  2—18  and  2—23)  for  a  solvent 
with  a  single  solvent  relaxation  time  rs  =  1.5  ps.  All  other  parameter  are  identical  to  the 
simulation  results  in  Figure  4  and  the  middle  panel  of  Figure  6.  At  long  times  ket(t)  is  well 
characterized  by  the  inverse  of  a  single  solvation  time,  with  a  limiting  rate  coefficient  ketm  = 
1/(1. 5  ps),  in  accord  with  simple  results  of  Section  3. 

The  overall  shape  of  ket(t),  as  time  increases,  is  qualitatively  similar  for  the  single 
exponential  and  bi— exponential  solvation  dynamics.  There  is  a  prompt  increase  in  ket(t)  at 
early  times  (barely  observable  in  Figure  6),  a  slower  decay  of  ket(t),  and  finally  a  limiting 
value.  The  early  time  increase  is  due  to  the  increase  in  the  magnitude  of  dF  (z)/dz  as  z 

o  1 

increases,  which  accelerates  the  diffusion.  The  intermediate  rate  coefficient  decrease  is 
apparently  due  to  the  strong  dependence  of  C(  (z)  on  z  as  shown  in  Figure  3. 


5.  The  Electron  Transfer  Kinetic  of  BA 

The  et  mechanism  of  BA  is  different  from  that  of  ADMA  in  two  important  ways. 
First,  there  are  two  CT  states  which  should  be  considered  (Figure  7).  Second,  the  adiabatic 
St  et  mechanism  in  polar  solvent  of  BA  involves  a  small  barrier  (AF^  ~  0.5  kcalmol'1) 
between  the  reactant  "well"  and  product  "wells"  (Figure  8),  while  for  ADMA  the  reaction  is 
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barrierless  in  Si. 

The  appropriate  equations  for  the  three  adiabatic  excited  BA  are  summarized  by  the 
following. 

I S i>  =  Cy>(z)|CT’>  +  C<  n(z)|LE>  +  C<^(z)|CT>.  (5-1) 

|S2>  =  C<  2>,(z)|CT>>  +  c<  2)(z)|LE>  +  C<|>(z)|CT>.  (5-2) 

|S3>  =  C<  3»(z)|CT’>  +  C«3>(z)|LE>  +  C<  3>(z)|CT>.  (5-3) 

For  BA  laser  excitation  at  400  nm  prepares  the  St  state  exclusively.  The 
time— dependence  of  p  (z,t)  has  been  simulated  as  shown  in  Figure  9.  The  upper  panel 
portrays  the  recently  published  GLE  simulation  of  p( z,t)  in  the  overdamped  limit  using  the 
following  solvation  parameters  for  C(t):  A!  =  0.46,  A2  =  0.54,  t\  —  0.43  ps,  and  r2  =  4.1  ps. 
p( z,t)  is  calculated  from  z(t)  trajectories  calculated  by  the  GLE.  The  lower  panel  of  Figure 
9  shows  a  GSE  simulation  of  p( z,t)  using  the  same  initial  conditions  and  solvation  parameter 
as  the  GLE  calculation.  The  evolution  of  p( z,t)  according  to  the  GSE  is  noticeably  slower 
(roughly  a  factor  of  2).  It  has  been  noticed  in  other  contexts  that  the  GSE  leads  to  slower 
(and  apparently  less  accurate)  dynamics.  In  the  present  case,  our  goal  is  to  explore  the 
relationships  between  ket(t),  A(t),  I(t,i/),  and  SG(t),  rather  than  explore  the  relationship 
between  ket(t)  and  the  solvation  dynamics,  so  we  emphasize  GbE  calculations  which  are 
computationally  less  expensive.  Note  that  the  small  barrier  between  the  reactant  and 
products  is  large  enough  to  maintain  a  reactant  peak  in  p  (z,t).  This  separation  of  p  (z,t) 
into  different  regions  coupled  with  the  S i/So  energy  gap  dependence  due  to  the  reaction 
barrier,  causes  the  simulated  spectra  to  be  approximately  the  sum  of  separate  "reactant" 
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and  "product"  bands,  see  Figure  10. 

The  presence  of  even  a  small  barrier,  as  in  the  case  of  BA,  is  apparently  sufficient  to 
dramatically  alter  the  spectral  manifestations  of  the  et  kinetics.  See,  in  particular,  in  Figure 
11  the  dynamics  of  the  survival  SQ(t)  ("exact",  see  eqns  2-18  and  2-23),  the  spectral  area 
A(t),  and  the  intensity  dynamics  I(t,413nm)  at  the  "blue  edge"  emission  wavelength.  All 
these  are  very  similar.  Thus,  either  experimental  observable  A(t)  or  I(t,^)  gives  a 
reasonably  accurate  measurement  of  S.(t).  For  example,  the  half— lifetime  of  S  (t)  in  Figure 
11  is  3.0  ps,  while  the  corresponding  values  for  A(t)  and  I(t,  4l3nm)  are  3.3  and  3.1  ps, 
respectively.  This  demonstrates  that  the  existence  of  even  a  small  barrier  can  simplify  the 
measurement  of  adiabatic  et  kinetics,  for  ket(t)  by  essentially  producing  a  "reactant"  region 
of  the  emission  spectrum.  This  is  analogous  to  the  large  barrier  limit  of  the  strongly 
adiabatic  regime,  which  resembles  transition— state— theory.  In  this  limit  the  survival  is 
defined  for  BA  as 

Stst  =  /  Cdz„s  (z,t),  (5-1) 

-Z  C 

where  zc  and  — zc  are  the  locations  of  the  dividing  points  separating  the  reactant  and 
products  wells.  In  an  earlier  paper  we  demonstrated  that  STST(t)  was  nearly 
indistinguishable  from  A(t)  and  I(t,413nm)  for  the  BA  simulations  using  the  strongly 
adiabatic  model. ^ 

The  fact  that  the  probability  distribution  has  separate  peaks  in  the  LE  and  CT 
regimes  (due  to  the  small  barriers  separating  the  LE  and  CT  minima  in  Fg^(z))  causes  the 

approximation  that  C.(  2(t)  ~  (C<  (t))2  to  be  invalid.  Consequently,  the  estimation  of 

ket(t)  from  A(t)  by  eqn  2-26  for  BA  is  in  error  and  completely  inappropriate. 
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From  the  standpoint  of  experimental  measurements  of  ket(t)  in  terms  of  the  strongly 
adiabatic  model,  it  is  essential  that  an  extensive  analysis  be  made  for  each  molecules  in 
order  to  classify  the  molecule  in  either  the  ADMA  class  (which  has  no  barrier  separating  the 
reactant  and  product  region)  or  the' BA  class  (which  possesses  a  barrier).  The  static 
absorption  and  fluorescence  spectroscopy  has  features  which  should  allow  for  a  rapid 
classification  for  molecules  that  fall  strongly  in  either  class.  ’  For  molecules  like  BA  the 
presence  of  the  barrier  ensures  that  S„(t)  can  be  quantified  directly  from  I(t,z/)  in  the  blue 
edge  of  the  fluorescence.  For  molecules  that  fall  in  the  ADMA  class,  A(t)  must  be  measured 
experimentally  in  order  to  estimate  S„(t)  (using  eqn  2—26)  or  alternatively  S  (t)  may  be 
calculated  "exactly"  using  eqns  2-18  and  2-23. 

5.  Conclusions  and  Summary 

The  recently  introduced  phenomenological  model  for  strongly  adiabatic  et  in  the 
excited  states  of  aromatic  electron  donor/acceptor  compounds  has  been  developed  in  more 
detail.  The  underlying  assumptions  in  the  model  have  been  delineated.  A  definition  of  a 
generalized  reactant  survival  probability  SG(t)  has  been  given.  An  expression  relating  SG(t) 
and  the  experimental  observable  A(t),  the  integrated  emission  intensity  has  been  derived. 
The  utility  of  these  approaches  have  been  applied  to  the  photodynamics  of  BA  and  ADMA, 
using  the  previously  estimated  free  energy  versus  solvent  coordinate  dependencies  of  these 
compounds  in  typical  organic  solvents.  The  general  applicability  of  the  strongly  adiabatic 
model  has  been  discussed. 
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Figure  Captions 

Figure  1.  Molecular  structures  of  bianthryl  (BA)  and 

4-{9-anthryl)-N,N’-dimethylaniline  (ADMA).  Donor  (D)  and  acceptor  (A)  signify  the 

electron  donor  and  acceptor  chromophores,  respectively.  Schematic  energy  levels  of  the 

three  relevant  states  in  the  non-polar  solvent:  ground  state  (DA  =  So),  locally  excited  state 
*  * 

(DA  =  LE)  which  results  from  the  r—r  excitation  of  the  anthracene  ring,  and  charge 
transfer  state  (D+— A-  =  CT)  which  is  an  intramolecular  ion— pair  state. 

Figure  2.  Theoretical  estimates  for  the  (left)  diabatic  and  (right)  adiabatic  free  energies 
of  ADMA  in  N,N’— dimethylformamide,  as  a  function  of  the  solvent  coordinate.  The  free 
energy  parameters  are:  ks  =  13.0  kcalmol'1,  AF°  =  10.5  kcalmol*1,  H  =  0.8  kcalmol'1. 

LLjLI 

See  ref.  3  and  4  for  details  of  the  simulation. 

Figure  3.  (upper  panel)  Theoretical  estimate  for  the  adiabatic  Si  surface  of  ADMA  in 
N,N ’-dimethylformamide.  (middle  panel)  The  square  of  the  coefficient  of  the  LE  character 
in  the  adiabatic  Si  state,  as  a  function  of  the  solvent  coordinate,  (lower  panel) 
Time-dependent  probability  distribution  function  on  the  Si  surface  of  ADMA  in 
N,N’— dimethylformamide  obtained  from  a  generalized  Smoluchowski  equation  (GSE)  and 
C(t)  (eqn  4—2)  with  Ai  =  0.55,  A2  =  0.45,  tx  =  0.75  ps,  and  7-2  =  2.5  ps  for  the  solvation 
dynamics.  See  ref.  3  for  the  details  of  the  GSE  simulation. 

Figure  4.  Simulation  of  the  time— resolved  emission  spectra  of  ADMA  in 
N,N’-dimethylformamide  obtained  from  the  probability  distribution  function  on  the 
adiabatic  S[  surface  shown  in  Figure  3.  The  ratio  of  the  transition  dipole  moment  (mct//*le) 
(eqns  2—28  and  2-29)  is  0.30  in  the  simulation. 


Figure  5.  (upper  panel)  The  integrated  emission  intensity  of  the  simulated 
time-resolved  spectra  (Figure  4)  obtained  from  eqn  2—26  for  ADMA  in 
N,N’-dimethylformamide  (DMF).  (middle  panel)  The  generalized  survival  (S  (t))  for  the 
adiabatic  et  reaction  on  the  Si  surface  of  ADMA  in  DMF  obtained  "exactly"  from  eqn  2—16. 
(lower  panel)  The  generalized  et  rate  coefficient  obtained  from  eqn  2-23  and  the  generalized 
survival  shown  above. 

Figure  6.  (upper  panel)  The  generalized  et  rate  coefficient  for  the  adiabatic  et  reaction 
on  the  Si  surface  of  ADMA  in  N,N’— dimethylformamide  obtained  "approximately"  from  the 
integrated  emission  intensity  and  eqn  2—33  with  the  approximation  of  eqns  2—30  and  2—31. 

A  two  exponential  function  (eqn  4—2)  with  Ai  =  0.55,  A2  =  0.45,  ri  =  0.75  ps,  and  r2  =  2.0 
ps  was  used  for  the  solvation  dynamics,  (middle  panel)  The  generalized  et  rate  coefficient 
obtained  "exactly"  from  the  generalized  survival  shown  in  Figure  5  and  eqn  2-23.  A  two 
exponential  function  with  the  same  parameters  as  above  was  used  for  the  solvation 
dynamics,  (lower  panel)  The  generalized  et  rate  coefficient  obtained  "exactly"  from  the  eqn 
2-23  and  the  probability  distribution  function  (eqn  2-14).  A  single  exponential  function 
with  rs  =  1.5  ps  was  used  for  the  solvation  dynamics. 

Figure  7.  Theoretical  estimates  for  the  (left)  diabatic  and  (right)  adiabatic  free  energies 
of  BA  in  propylene  carbonate  as  a  function  of  the  solvent  coordinate.  The  free  energy 
parameters  are:  ks  =  21.0  kcalmol'1,  AF°  =  5.0  kcalmol'1,  and  Hrp  rT  =  1.0  kcalmol'1.  See 
ref.  3  and  4  for  details  of  the  simulation. 

Figure  8.  (upper  panel)  Theoretical  estimate  for  the  adiabatic  S!  surface  of  BA  in 
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propylene  carbonate  as  a  function  of  the  solvent  coordinate,  (lower  panel)  The  square  of 
the  coefficient  of  the  LE  character  in  the  adiabatic  St  state.  The  free  energy  parameters  for 
both  the  panels  are  the  same  as  those  in  Figure  7. 

Figure  9.  (upper  panel)  Time-dependent  probability  distribution  function  p  (z,t)  on 

5 1 

the  Si  surface  of  BA  in  propylene  carbonate  obtained  from  a  generalized  Langevin  equation 

(GLE)  and  C(t)  (eqn  4—2)  with  Aj  =  0.46,  A  2  =  0.54,  ri  =  0.43  ps,  and  =  4.1  ps  for  the 

solvation  dynamics,  (lower  panel)  p  (z,t)  for  the  same  chemical  obtained  from  a 

**  1 

generalized  Smoluchowski  equation  (GSE)  and  the  same  functional  form  as  above  for  the 
solvation  dynamics.  See  refs.  4  and  3  for  details  of  the  GLE  and  GSE  simulations, 
respectively. 

Figure  10.  Simulation  of  the  time— resolved  emission  spectra  of  BA  in  propylene 

carbonate  obtained  from  the  distribution  function  on  the  St  surface  shown  in  the  lower  panel 
of  Figure  8.  A  single  exponential  function  of  C(t)  (eqn  4-2)  with  rs  =  2.4  ps  was  used  for 
the  solvation  dynamics.  The  ratio  of  the  transition  dipole  moments  (p  /p  )  (eqns  2—28 
and  2—29)  is  0.32  in  the  simulation. 

Figure  11.  (upper  panel)  The  integrated  emission  intensity  of  the  simulated 
time-resolved  spectra  (Figure  10)  obtained  from  eqn  2—26  for  BA  in  propylene  carbonate 
(PC),  (middle  panel)  The  generalized  survival  (SG(t))  for  the  adiabatic  et  reaction  on  the 
Si  surface  of  ADM  A  in  PC  obtained  "exactly"  from  eqn  2—18.  (lower  panel)  The  transient 
decay  dynamics  of  the  simulated  emission  spectra  (Figure  10)  at  413  am  ("the  blue  edge"). 
For  all  the  panels  a  single  exponential  function  of  C(t)  (eqn  4—2)  with  rs  =  2.4  ps  was  used 
to  represent  the  solvation  dynamics. 
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